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Abstract. A convection-driven MHD dynamo in a rotating spherical shell, with clearly 
defined structural elements in the flow and magnetic field, is simulated numerically. Such 
dynamos can be called deterministic, in contrast to those explicitly dependent on the 
assumed properties of turbulence. The cases most interesting from the standpoint of 
studying the nature of stellar magnetism demonstrate the following features. On a global 
scale, the convective flows can maintain a "general" magnetic field with a sign-alternating 
dipolar component. Local (in many cases, bipolar) magnetic structures are associated with 
convection cells. Disintegrating local structures change into background fields, which drift 
toward the poles. From time to time, reversals of the magnetic fields in the polar regions 
occur, as "new" background fields expel the "old" fields. 



Introduction 

The original motivation of this study was suggested by the 
problems of solar physics. Although the results obtained at 
this stage cannot be interpreted to represent the specifically 
solar dynamo process, we should give a brief exposition of 
the previous investigation that led us to the formulation of 
the problem considered here. 

Observations of the solar magnetic fields reveal a bewil- 
dering variety of structures and activities. It is remarkable 
that solar processes vary in scale from sizes comparable to 
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the solar radius to the limit of present resolution and in 
duration from tens of years to minutes (see, for instance, 
Schrijver and Zwaan [2000]). 

Mean-field electrodynamics [Mojfatt, 1978; Krause and 
Rddler, 1980] has clarified many issues concerning the gen- 
eration of the global magnetic fields of cosmic bodies. How- 
ever, such problems as the formation of local magnetic fields 
and their relationship to the global fields fall completely be- 
yond the scope of mean-field theories. Meanwhile, the phe- 
nomenon of solar and stellar magnetism could be adequately 
understood only if the dynamics of the interplay between 
structures in the velocity field and magnetic field is compre- 
hensively studied over a wide range of spatial scales. Dy- 
namo models that are aimed at a unified description of the 
global and local processes and that deal with local, instan- 
taneous quantities rather than averaged ones can naturally 
be referred to as "deterministic" models. They describe the 
structural elements present in the flow and in the magnetic 
field instead of considering the averaged parameters of the 
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turbulent flow (in particular, the statistical predominance of 
one sign of the velocity-field helicity or another). 

The idea that convection cells in the solar subphoto- 
spheric zone could be a connecting link between global and 
local magnetic fields traces back to the mid-1960s. Tver- 
skoy [1966] represented the convection cell by a toroidal 
eddy and demonstrated, in the framework of a kinematic 
approach, that such a model convection cell can amplify the 
magnetic field and produce characteristic bipolar magnetic 
configurations. This approach was also used by Getling and 
Tverskoy [1971] to construct a kinematic model of the global 
dynamo in which toroidal eddies distributed over a spheri- 
cal shell, acting jointly with the differential rotation of the 
shell, maintain a sign-alternating global magnetic field. If 
a poloidal magnetic field is present, the differential rotation 
produces a toroidal component of the global magnetic field. 
If the local magnetic configuration produced by an eddy in- 
teracting with the large-scale toroidal field is rotated through 
some angle about the axis of the eddy, this configuration 
contributes to the regeneration of the poloidal component 
of the global magnetic field. Thus, a cell locally interacting 
with the magnetic field serves in this model as a building 
block of the global dynamo, and the latter can in this case 
be called the "cellular" dynamo. The rotation of the local 
magnetic-field pattern can be expected if the system rotates 
as a whole and the flow is affected by the Coriolis force. 

In recent years, after the advent of suitable comput- 
ing facilities, some steps have been made to verify these 
ideas by means of numerical simulation. Getling [2001] and 
Getlmg and Ovchmnikov [2002] obtained numerical solutions 
to the three-dimensional nonlinear problem of magnetocon- 
vection in a plane horizontal layer of incompressible fluid, 
heated from below, and found that hexagonal convection 
cells interacting with a weak initial ( "seed" ) , horizontal mag- 
netic fleld can produce various structures of the strongly 
amplified magnetic field, with a predominant bipolar com- 
ponent. Dobler and Getling [2004] extended this numerical 
analysis to compressible fluids and obtained similar results. 

Modern computing resources make it possible to approach 
the development of numerical cellular-dynamo models that 
could provide a parallel description of both the global and lo- 
cal magnetic flelds. However, even today, numerical schemes 
can hardly be used to simulate flows and magnetic flelds 
over scale ranges covering two or more orders of magnitude. 
Only the largest convection cells — of sizes comparable to 
the depth of the convection zone (such as solar "giant" cells, 
for which little observational evidence exists) — can be sim- 
ulated in the framework of global models. If we assume that 
the principal features of the process should be similar for 
convection on different scales, such global models would help 
us to verify our qualitative notion and provide guidelines for 
the elaboration of a more detailed description. 

Here, we use numerical simulations to investigate the 
properties of cellular dynamos in rotating spherical shells, 
which could operate in stars under certain conditions. Al- 
though there are reasons to believe that the solar dynamo is 
also of the cellular type, we do not directly associate the 
presently obtained results with solar processes, since the 
computed patterns of magnetic-filed evolution bear only lim- 
ited similarity to the pattern observed on the Sun. We 




Figure 1. Geometrical configuration of the problem. The 
two spherical surfaces with radii n and ro are dark-shaded. 
A part of the outer surface is removed to expose the inte- 
rior of the shell (light-shaded) where the conducting fiuid is 
confined. 



merely note some remarkable features of the cellular dy- 
namos, which may be of interest from the standpoint of 
stellar magnetohydrodynamics. 



Formulation of the problem and numerical 
technique 



In order to model the process of magnetic-field generation 
in a stellar convection zone, we consider a spherical shell of 
thickness d — ro — ri (where ro and ri are the outer and 
inner radii of the shell), full of electrically conducting fiuid 
and rotating with a constant angular velocity ft about a 
fixed axis e^, as shown in Figure 

We follow the standard formulation used in earlier work 
by Tilgner and Busse [1997], Busse et al. [1998], Grote et 
al. [1999, 2000], and Simitev and Busse [2002, 2005], but 
we assume a more general form of the static temperature 
distribution, 



Ts = To ~ -d r + —-, 
2 d r 



(1) 



Pi = 



■qdAT 
(1^ 



where the radial coordinate r is measured in units of d, x 
is the thermal diffusivity, Cp is the specific heat at constant 
pressure, q is the mass density of uniformly distributed heat 
sources, rj = ri/ro is the inner-to-outer radius ratio of the 
shell, and To is a constant. The quantity AT is related 
to the difference between the constant temperatures of the 
inner and outer spherical boundaries, Ti and To, as 



AT = Ti- To 



2 1 — 



(2) 
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and reduces to Ti — To in the case of q — (also dealt with 
in some simulations). The shell is self-gravitating, and the 
gravitational acceleration averaged over a spherical surface 
r — const can be written as g = —'yd ■ r, where r is the 
position vector with respect to the centre of the sphere; as 
specified above, its length r is measured in units of d. In ad- 
dition to d, the time d^ /u, the temperature v'^/jad'^ (where 
a is the volumetric coefficient of thermal expansion), and 
the magnetic induction u{jj,pY^^ / d are used as scales for the 
dimensionless description of the problem; here, v denotes 
the kinematic viscosity of the fiuid, p is its density, and /i is 
its magnetic permeability (we set p= 1). 

We use the Boussinesq approximation in that we assume 
p to be constant except in the gravity term, where, in ad- 
dition to the standard linear dependence p{T) [according to 
which p~^ (dp/dT) = —a = const], we introduce a small 
quadratic term in most cases. Once a cellular pattern has 
developed, the presence of this term and of the volumetric 
heat sources should not radically modify the properties of 
the dynamo; however, both these factors favour the devel- 
opment of polygonal convection cells [Busse, 2004] similar 
to the cells observed on the Sun, rather than meridionally 
stretched, banana-like convection rolls. Without these es- 
sential modifications, polygonal cells could only be obtained 
at much smaller rotational velocities; in this case, the pro- 
cess would develop very slowly, and the computations would 
be extremely time-consuming. 

Thus, the equations of motion for the velocity vector u, 
the heat equation for the deviation O from the static tem- 
perature distribution, and the equation of induction for the 
magnetic field B are 



V-u = 0, 

{dt + U- V)U : 



-Vtt -I- ru X -f (e -f- e 6 )r 
-f V^u+ (V X B) X B, 



P{dt + u-V)Q = V^e + [n + R. 



V-B = 0, 

atB = V X (u X B) 



-P^'V'B, 



(3a) 

(3b) 

u, 

(3c) 
(3d) 
(3e) 



where vr is an effective pressure. 

Six nondimensional physical parameters of the problem 
appear in our formulation. The Rayleigh numbers measure 
the energy input into the system. 



Rc 



(4) 



and are associated with the internally distributed heat sour- 
ces q and the externally specified temperature difference Ti — 
To [see equation Q], respectively. The Coriolis number r, 
the Prandtl number P, and the magnetic Prandtl number 
Pm describe ratios between various time scales in the system. 



2nd' 



(5) 

I' X J^m 

{i/m is the magnetic viscosity, or magnetic diffusivity) . Fi- 
nally, e is the small constant that specifies the magnitude of 
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Figure 2. Static profiles of the temperature and tempera- 
ture gradient in the case of internal heating. 



the quadratic term in the temperature dependence of density 
[see equation Il3bll 1. 

Since the velocity field u and the magnetic induction B 
are solenoidal vector fields, the general representation in 
terms of poloidal and toroidal components can be used. 



u = V X (Vw X r) -f V™ X r, 
B = V X (V/i X r) + Vg X r. 



(6a) 
(6b) 



By multiplying the (curl)^ and the curl of the Navier-Stokes 
equation Il3bt in the rotating system by r, we obtain two 
equations for v and w, 

[(V' - dt)£2 + rd^]V\ + tQw - £2(6 -f- e 8^) 

= -r • V X [V X (u • Vu - B • VB)] 

(7a) 

[(V^ - dt)C2 + Td^]w - tQv = r • V X (u ■ Vu - B ■ VB), 

(7b) 

where ip denotes the azimuthal angle ("longitude") in the 
spherical system of coordinates r, 0, ip, and the operators £2 
and Q are defined by 

£2 = -rV + 9,(r'a,) 

Q = rcosOV^ — (£2 -I- rdr){cosOdr — smQde). 
The heat equation Il3cll can be rewritten in the form 

V^e + [n + i?er7r-^(l - r?)-^] £2^; = P(dt + u ■ V)e. (8) 

Equations for h and g can be obtained multiplying the equa- 
tion of magnetic induction I3el and its curl by r, 

V^£2/i = Pm[at£2/i-r ■ V X (u X B)], (9a) 
V^£2g = Pm[dtC2g - r • V x ( x(u x B))]. (9b) 

We assume stress-free boundaries with fixed tempera- 
tures, 

V — drrV = dr{w/r) — Q = Q aX r = r\ and r — To- (10) 
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For the magnetic field, we use electrically insulating bound- 
aries such that the poloidal function h must be matched to 
the function h^-^^ that describes the potential fields outside 
the fluid shell 

g = h- = dr{h - ft''^') = at r = n and r- = r^. (11) 

The numerical integration of equations ffl- llllll proceeds 
with a pseudospectral method developed by Tilgner and 
Busse [1997] and Tilgner [1998] , which is based on an expan- 
sion of all dependent variables in spherical harmonics for the 
9 and Lp dependences; in particular, for the magnetic scalars, 

^ oo I 

9=^Y. E Gr(r,i)PrWexp{im^}, (12a) 

- oc I 

'^=-E E Hrir,t)Prie)eMimp} (12b) 

1=0 m=-l 

(with truncating the series at an appropriate maximum /), 
where P/" denotes the associated Legendre functions. For 
the r dependences, truncated expansions in Chebyshev poly- 
nomials are used. The equations are time-stepped by 
treating all nonlinear terms explicitly with a second-order 
Adams-Bashforth scheme whereas all linear terms are in- 
cluded in an implicit Crank-Nicolson step. 

For the computations to be reported here, a minimum 
of 33 collocation points in the radial direction and spheri- 
cal harmonics up to the order 96 have been used. In addi- 
tion to the geometric parameter 77 and the above-mentioned 
physical parameters, we specified a computational parame- 
ter, viz., the fundamental (lowest nonzero) azimuthal num- 
ber mo. Thus, only the following azimuthal harmonics were 
really considered: 

^ g±i™0'i= g±2imov g±3imoip 




Figure 3. Contours of the radial velocity component on 
the sphere r = ri -I- 0.5 a± t = 98.73 (top) and of the radial 
component of the magnetic field on the sphere r = ro aX t = 
98.73 (middle) and 101.73 (bottom) in the case of internal 
heating at 77 = 0.6, Ri = 3000, i?c = -6000, r = 10, P = 1, 
Pm = 30, and mo = 5. Solid (red) curves: positive values; 
dotted (green) curves: zero values; dashed curves (blue): 
negative values. 



In other words, we imposed an mo-fold symmetry in the tp 
direction. If mo 7^ 1, this reduces the computation time. 



Results 
Internal heating 

For the cases of internal heating, we varied Pm and as- 
sumed 77 = 0.6, Pi = 3000, Pc = -6000, r = 10, P = 1, 
and mo = 5. As can be seen from checking computations 
with mo = 1 (not presented here), removing the artifi- 
cially imposed fivefold azimuthal symmetry does not sub- 
stantially modify the character of the convection pattern. 
The quadratic term was present in the temperature depen- 
dence of density, with a control parameter of e = 0.005. 

The distributions of the temperature Ta,{r) and its gra- 
dient dTs/dr for the corresponding static-equilibrium state 
are shown in Figure|5| Obviously, the outer part of the shell 
is convectively unstable and the inner part is stable. 




Figure 4. Contours of the axisymmetric azimuthal velocity 
(left-hand section) and streamlines of the meridional circula- 
tion, or contours of the stream function of the axisymmetric 
meridional flow (right-hand section) in the same case with 
internal heating (Pm = 30) at t = 98.73. The curves have 
the same meaning as in Figure 3. 
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Figure 5. Variation in the amplitude of the dipolar har- 
monic of the poloidal magnetic field, Hl{r, t), at r = ri + 0.5 
in the same case with internal heating (Pm = 30). 



The case of Pm = 30. At this Pm value, the compu- 
tations covered a time interval of about 100 in units of the 
time of thermal diffusion across the shell. Over most part 
of this period, a very stable pattern of convection cells with 
a dodecahedral symmetry can be observed (Figure|3^ top). 
These cells have a normal appearance typical of cellular con- 
vection, without substantial distortions due to the rotation 
of the shell. The entire pattern drifts in the retrograde direc- 
tion, in agreement with theoretical predictions [Busse, 2004]. 

The axisymmetric component of the azimuthal velocity 
(Figure [IJ in a well-established flow pattern is nearly sym- 
metric with respect to the equatorial plane. Specifically, a 
prograde rotation of the equatorial zone (in the frame of 
reference rotating together with the entire body) is present 
along with a retrograde rotation of the midlatitudes, and 
pairs of "secondary" prograde- and retrograde-rotation zones 
can also be noted in the polar regions. In a nonrotating 
frame of reference, the equatorial zone rotates more rapidly 
and the midlatitudinal zones more slowly than the shell as 
a whole does. Three pairs of meridional-circulation vortices 
fill the entire meridional section of the shell, from one pole 
to another. 

The pattern of magnetic field is less regular than the pat- 
tern of fiow (Figure |3] middle and bottom). Some remark- 
able features or the simulated dynamo process can be sum- 
marized as follows. 

First, local magnetic structures associated with convec- 
tion cells emerge repeatedly as compact magnetic regions 
(see Fig. I^J. In their subsequent evolution, these regions 
change their configuration and finally dissipate into much 
weaker remnant fields. 

Second, the dipolar component of the global magnetic 
field exhibits polarity reversals [see Figure |S] for a graph of 
the amplitude of the dipole component, _ffi(t)]^. The back- 
ground fields — remnants of the decaying local magnetic 
structures — drift toward the poles and "expel' the "old" 
background fields present in the polar regions. As a result, 
the old magnetic polarity is replaced with the new one due to 
the poleward drift of the latter. The polarity reversals of the 
global magnetic field can also be seen from the variation in 
the amplitude of the dipolar harmonic of the poloidal field. 



^ The colour figures can be viewed in the online version of the 
paper. 

^Wherever the r variable as an argument of is omitted, we 
mean //^(r-i -1-0.5, t) 



Hlit) (FigureEJ. The two maps of the magnetic field shown 
in Figure |3| correspond to two situations in which the global 
magnetic dipole has opposite orientations (the polarity re- 
versal between these two times corresponds to the rightmost 
intersection of the curve in Figure|^with the horizontal zero 
line). 

Third, an interesting intermittent behaviour is exhibited 
by the magnetic energy of the system. Let us compare the 
full energy and two particular fractions of the energy associ- 
ated with the magnetic-field component that has a dipolar- 
type symmetry (i.e., is antisymmetric with respect to the 
equatorial plane). Specifically, we are interested in the be- 
haviour of the energy of the axisymmetric and the nonax- 
isymmetric part of this component. The axisymmetric part 
is represented by the spherical harmonics with I odd and 
m = Q [see (1121 1. and the nonaxisymmetric part by other 
harmonics with / 4- m odd. As can be seen from Figure |S| 
(in which the total energy and its particular fractions are 
divided by the volume of the shell), the main peaks in the 
graph of the total energy are alternately associated with in- 
creases in the energies of the axisymmetric and the nonax- 
isymmetric part of the component with a dipolar symmetry. 
In particular, the peak located near t = 42 is fed by the 
symmetric field; near t = 55, by the asymmetric field; near 
i = 68, by both but with some predominance of the sym- 
metric part; and near t = 78, again by the asymmetric part. 

The same case of Pm = 30 but with special initial 
conditions. In our attempts to find conditions for the re- 
alization of magnetic-field dynamics similar to the generally 
imagined pattern of a hypothetical dynamo process with 
differential rotation as its essential part (known since the 
qualitative model suggested by Babcock [1961] and Leighton 
[1964, 1969]), we made an additional computational run. 
We specified all parameters to be the same as in the above- 
described case. However, the initial conditions were chosen 
in such a way that, initially, the system would more likely 
find itself within the attraction basin of the expected dy- 
namo regime in state space. To this end, we superposed the 
nonaxisymmetric components of the velocity field and mag- 
netic field obtained in the above-described run onto a pattern 
of differential rotation with the equatorial belt accelerated 
compared to higher latitudes and with an appropriate dis- 
tribution of the axisymmetric azimuthal magnetic field that 
has different signs on the two sides of the equator. 

The computations have demonstrated that the system 
nevertheless approaches virtually the same regime that was 
observed without using such special initial conditions. Thus, 
a closer similarity between the numerical solution and the 
properties of the hypothetical dynamos of the Babcock- 
Leighton type does not seem to be achievable in the frame- 
work of this very simple model. 

The case of Pm = 200. As the magnetic Prandtl num- 
ber Pm is varied (under otherwise fixed conditions), the con- 
vection pattern varies little over a fairly wide Pm range. 
However, the greater this parameter, the higher the mean 
strength of the magnetic field (and, accordingly, the total 
magnetic energy). The kinetic energy of convection de- 
creases with the increase of Pm and convection becomes more 
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Figure 6. Density of magnetic energy and its particular fractions in the same case with internal heating (Pm = 30). Solid 
(heavy black) curve and right-hand vertical scale: total density of the magnetic energy; dashed (red) curve and left-hand 
vertical scale: energy density associated with the axisymmetric part of the magnetic-field component that has a dipole 
symmetry; dotted (blue) curve and left-hand scale: energy density associated with the nonaxisymmetric part of the same 
magnetic-field component. 




Figure 7. Contours of the radial component of the velocity on the sphere r = ri 0.5 (left) and of the magnetic field on 
the sphere r = ro (right) a,t t — 200 (top) and t = 327.2 (bottom) in the case of internal heating at = 200, while the 
other parameters are as in the preceding case: rj = 0.6, Ri = 3000, Rc = —6000, r = 10, P = 1, and mo = 5. The curves 
have the same meaning as in Figure 3. 



sensitive to time variations in the magnetic field. The in- 
crease of Pm is also manifest in the fact that local magnetic 
fields become more patchy and less ordered. Individual ar- 
eas filled with the magnetic field of a given sign are smaller 
in size and more numerous, and bipolar structures are not 
so well pronounced (see Figure[7| which refers to Pm ~ 200). 
As in the case of Pm = 30, we can observe the penetration 
of background fields into the polar regions and sign reversals 
of the polar background fields. 

Our computation for Pm = 200 cover a time interval al- 
most five times as long as for Pm = 30 (Figure |HJ. The two 
velocity maps and two magnetic-field maps shown in Fig- 
ure|7]nearly correspond to the times of one negative and one 
positive extremum of the amplitude Hi(t) (see Figure |HJ. It 
is remarkable that the polar background fields have different 
polarities at these two times. At t = 200.0, the background 
magnetic field is negative in the "northern" and positive 



in the "southern" polar region; an opposite situation takes 
place at t = 327.2. The H^{t) curve demonstrates numerous 
sign reversals, although fine details of this dependence only 
reflect the irregular, fluctuational aspect of the process. It 
is nevertheless clear that, even if we smooth this curve, it 
will exhibit quite pronounced cyclic, although nonperiodic, 
polarity reversals of the dipolar component of the "general" 
magnetic fleld. 

It should be noted that the distribution of the axisym- 
metric component of the azimuthal velocity (the pattern of 
differential rotation, not shown here) in this case is much 
more complex and variable than at Pm ~ 30. This effect 
also can be due to the stronger influence of the magnetic 
field on the fluid motion. 

Internal heating without a ©'^ term (nonmagnetic 
case). To form an idea of the role played by the quadratic 
term in the p{T) dependence, we computed a purely hydro- 
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Figure 8. Variation in amplitude of the dipolar harmonic of the poloidal magnetic field, H^{r,t), at r = n + 0.5 for the 
same case with internal heating as illustrated by Figure |7|(Pm = 200). 



dynamic (with B = 0) scenario under conditions that dif- 
fered from the conditions of the above-described simulation 
by the absence of the term (e = 0) and by the Coriolis 
number (r = 1); in addition, we assumed mo = 1 in this 
case. 

The principal result of these computations is the finding 
that, in the absence of the quadratic term, convection unaf- 
fected by the magnetic field forms patterns of well-localized, 
three-dimensional cells, which typically appear as shown in 
Figure |5| A downwelling is observed in the centre of each 
cell, in contrast to the above-described cases, in which cen- 
tral upwellings typically developed. The issue of the direc- 
tion of circulation in a convection cell is a fairly subtle matter 
(see, e.g., Getling [1999] for a survey of some situations re- 
lated to convection in horizontal layers), so that agreement 
or disagreement between our model and any really observed 
pattern can in no way be an indication for the factors respon- 
sible for the observed direction of convective motions. Our 
primary interest in the cases where an upwelling is present 
in the central part of a cell is merely dictated by our inten- 
tion of constructing a dynamo model reproducing the solar 
phenomena as closely as possible. 




Figure 9. Contours of radial velocity on the surface r — 
n + 0.5 in the nonmagnetic case with internal heating, r = 1; 
no quadratic term is present in the p(T) dependence (e = 0). 
The other parameters are i] — 0.6, Ri — 3000, -Rc ~ —6000, 
P = 1 (as in the preceding cases), and mo — 1. The curves 
have the same meaning as in Figure 3. 



Heating "from outside" (through the inner surface) 

In addition, we undertook a search for regimes in which 
convection preserves its "three-dimensional" structure in the 
absence of internal heat sources. In other words, some com- 
putations were done at -Ri = 0. The quadratic term was also 
missing from the temperature dependence of density (e = 0) 
in these runs. 

Note that, in the limiting case of a nonrotating shell, con- 
vection cells are not stretched in any direction. Therefore, a 
cellular pattern of convective motion can obviously be main- 
tained even without such favourable factors as a specific form 
of stratification and a quadratic term, but at smaller rota- 
tional velocities. 

In shells without internal heat sources {q — 0), convec- 
tion regimes similar to those observed at g 7^ should be 
expected, under otherwise identical conditions, at smaller rj. 
This is because a stratification similar to that shown in Fig- 
ure |5| confines the development of convective motion to the 
outer part of the shell. 




Figure 10. Map of the radial velocity component on the 
sphere r — ri + 0.5 in the nonmagnetic case with heating 
through the inner boundary at 77 = 0.9, P — 1, t — 0.1, 
-Rc ~ 5000, and mo = 2. The contours merge into hatching: 
moderately light (red) for positive values, light (green) for 
zero values, and dark (blue) for negative values. 
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We illustrate here the case of g = only by a tentative 
computational run for convection without a magnetic field, 
at = 0.9, P = 1, r = 0.1, = 5000, and m-o = 2 (see 
Figure nUl for a velocity field typical of this case). Although 
the convection pattern is complex in this case, a tendency 
toward the formation of meridionally elongated cells can nev- 
ertheless be noted. 

Computations with the magnetic field included (e.g., for 
ri = 0.8, P = 1, i?c = 5000, T = 0.1, = 5) demonstrate 
the development of magnetic features with a very small spa- 
tial scale, close to the resolution limit of the computational 
scheme. This is a signature for an insufficient spatial reso- 
lution, so that the simulation results are not quite reliable. 

The qualitative aspects of the results suggest that the 
Coriolis number proves again to be insufficiently small for 
the stability of "three-dimensional" convection cells, and the 
cells ultimately become substantially stretched, although not 
in a strictly meridional direction. 

On the whole, the last two scenarios of flow and magnetic- 
fleld evolution indicate that regimes of "cellular" dynamo in 
a shell without internal heating and without a quadratic 
term in the p{T) dependence should be sought in the range 
of smaller Q (and r). 



Conclusion 

We have constructed relatively simple numerical models 
that describe a self- sustained process of generation of inter- 
acting global and local magnetic fields. As in most hypo- 
thetical astrophysical dynamos, the generation is driven by 
thermal convection in combination with differential rotation. 
We did not introduce any kinematic elements in our model, 
so that the entire velocity field appeared as the solution of 
the full system of MHD equations. 

The most remarkable features revealed in the computed 
dynamo regimes can be summarized as follows. The pro- 
cess of magnetic-field generation is cyclic, although rather 
irregular. It includes the repeated generation of local, in 
many cases bipolar, magnetic structures. These structures 
dissipate giving rise to chaotic background fields. They may 
drift in the poleward direction, replacing the already exist- 
ing, "old" background fields. In some cases, a correspon- 
dence can be noted between such polarity reversals in the 
polar regions and sign reversals of the axisymmetric bipolar 
component of the global magnetic field. 

One of the computed scenarios demonstrates a remarkable 
intermittency in the behaviour of some fractions of magnetic 
energy: the axisymmetric and the nonaxisymmetric part of 
the magnetic-field component with a dipolar symmetry al- 
ternate in making larger contributions to the total-energy 
peaks. 

Mean-field dynamo models, which have been most popu- 
lar in astrophysics over a few past decades, attribute the gen- 
eration of the global magnetic fields of stars to the a effect - 
the statistical predominance of one sign of the velocity-field 
helicity over another. It is quite plausible that the a effect, 
in one form or another, is a fairly general property of vari- 
ous velocity fields capable of maintaining undamped regular 



magnetic fields. However, this property must not necessar- 
ily be associated with turbulent motion. In particular, the 
model velocity field in the toroidal eddies used by Getling 
and Tverskoy [1971] to construct a global dynamo included 
an azimuthal (with respect to the axis of the eddy) veloc- 
ity component, so that the trajectories of the fluid particles 
were spirals deformed in a certain way. A similar property 
may also be inherent in the convective flow that develops 
in our model, although checking this possibility requires a 
special investigation. 

At this stage, the "deterministic" cellular dynamo de- 
scribed here is oversimplified to be regarded as a model of 
the solar or any other specific stellar dynamo. However, it 
demonstrates that the dynamics of the well-structured lo- 
cal magnetic fields and of the global magnetic field may be 
ingredients of one complex process. 
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